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ABSTRACT 

The Multi-Epoch Nearby Cluster Survey (MENeaCS) has discovered twenty-three cluster Type la 
supernovae (SNe) in the 58 X-ray selected galaxy clusters (0.05< z <0.15) surveyed. Four of our 
SN la events have no host galaxy on close inspection, and are likely intracluster SNe. Although one of 
the candidates, Abell399_3_14_0, appears to be associated in projection with the outskirts of a nearby 
red sequence galaxy, its velocity offset of ~ 1000 km s _1 indicates that it is unbound and therefore 
an intracluster SN. Another of our candidates, Abell85_6_08_0, has a spectrum consistent with a SN 
1991bg-like object, suggesting that at least some portion of intracluster stars belong to an old stellar 
population. Deep image stacks at the location of the candidate intracluster SNe put upper limits 
on the luminosities of faint hosts, with M r > —13.0 mag and M g > —12.5 mag in all cases. For 
such limits, the fraction of the cluster luminosity in faint dwarfs below our detection limit is < 0.1%, 
assuming a standard cluster luminosity function. All four events occurred within ^600 kpc of the 
cluster center (projected), as defined by the position of the brightest cluster galaxy, and are more 
centrally concentrated than the cluster SN la population as a whole. After accounting for several 
observational biases that make intracluster SNe easier to discover and spectroscopically confirm, we 
calculate an intracluster stellar mass fraction of 0.16^; Jg (68% confidence limit) for all objects within 
i?200- If we assume that the intracluster stellar population is exclusively old, and the cluster galaxies 
themselves have a mix of stellar ages, we derive an upper limit on the intracluster stellar mass fraction 
of < 0.47 (84% one-sided confidence limit). When combined with the intragroup SNe results of McGee 
& Balogh, we confirm the declining intracluster stellar mass fraction as a function of halo mass reported 
by Gonzalez and collaborators. 



1. INTRODUCTION 
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Galaxy clusters contain a population of stars gravita- 
tionally unbound to individual galaxies, yet still bound 
to the clusters' overall gravitational potential, created 
by the stripping of stars from galaxies during interac- 
tions and mergers. Much work has been done to char- 
acterize the distribution, stellar mass and stellar pop- 
ulation of the intracluster light ( ICL) through a vari- 
ety of observational tracers (e.g. [Gal- Yam et all 120031: 
l Arnaboldi et al.l[2004l: iZibetti et alll2005llGonzalez et al l 
20051 iKrick fc Bernstein! 120071 : I Williams et all 120071 : 
iRudick et all 120101 among others) because it is an 
important remnant of the hierarchical accretion his- 
tory of clusters and the environmental processes which 
can transform cluster galaxies. Numerical work has 
sought to reproduce the observed results, to under- 
stand the processes that form the ICL, and how those 
processes in turn a ffect cluster galaxies and the clus- 
ter a s a whole (e.g. iMurante et al1l2004t iWillman et all 



1 20041: iSommer-Larsen et al l 120051 ; iPurcell et all 120081: 
iPuchwein et al.ll2010l iDolag et all 20101 among others) . 

One independent and complementary technique for 
probing the ICL a re hq stless, IC SNe. For instance, 
iSivanandam et all (|2009( ) have shown that, under rea- 
sonable assumptions, IC SNe can generate ^30% of the 
iron observed in the IC medium out to rsoo, although 
this is subject to the still-large uncertainties in the SN 
la rate in galaxy clusters (especially in the relative rates 
of hostless and hosted SNe) and the evolution of the 
ICL with respect to the cluster as a whole. More rel- 
evant for the current work, the relative numbers of host- 
less and hosted SNe la have been used as a measure of 
the IC stellar mass fraction (e.g. iGal-Yam et all 120031 : 
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iMcGee fc Baloghll2010l) . Additionally, hostless SNe can 
serve as a tracer of the extent and distribution of the 
ICL - even the most diffuse and low surface brightness 
elements of the ICL can be revealed (albeit with relative 
rarity) if a SN explodes in its midst. As IC SNe observa- 
tions become more common and better studied, the rel- 
ative mix of core collapse SNe and SNe la can constrain 
the age of the underlying stellar population. Indeed, 
even the relative mix of the known sub-classes of SNe 
la - where SN 1991T-like events (and others with slow- 
declining light c urves) are associated with young stellar 
populations (e.g. lHamuv et al.lll995[ ). while SN 1991bg- 
like events (and those with fast-d eclining light curves) 
are associated with old ones (e.g. IHowelll [20011 - can 
serve as a tracer of the underlying stellar population. 

Intracluster SN observations are still rare. A fi rst pos- 
sible IC SN la was discovered by iSmithl (|1981h . in an 
apparent stream connectin g NGC 4406 and NG C 4374 
in the Virgo cluster. iGal-Yam et al.l (|2003t ) found 
two IC SNe out of a total sample of seven in the 
WOQTS survey, a c l uster SN search center ed at z ~ 0.1 
(iSharon et al.l 120071: lG~al-Yam et al.l l2008f ). More re- 
centlv. IMcGee fc Baloghl (|2010| ) mined the Sloan Digital 
Sky Survey (SDSS) SN Survey for SNe la associated with 
low redshift galaxy groups. Of their 59 group SNe la, 19 
have no detectable host galaxy, and after making correc- 
tions for the group luminosity fu nction and assuming a 
two-component SNe la rate (e.g. iMannucci et al.1 [20051 : 
iScannapieco fc Bildstenl [2001 they found that 47±^% 
of the stellar mass in thei r groups was in the form of 
diffuse light. Additionally, iDildav et all (|2010al ) discov- 
ered up to three hostless, IC SNe la in their study of the 
SNe la rate in galaxy clusters with the SDSS-II SN Sur- 
vey. Several other possible IC SNe have been discovered 
over the years (e . g. iGermanv 1998; Rciss & S abinelll998l: 
Sand et al.l [20081: IGal-Yam fc Simonll2008t ISharon et"aLl 



2010: IBarbarv et al.ll2010i r As outlined in the previous 



paragraph, hostless SNe are a unique probe of the diffuse 
universe. Much more can be done, and there is a need 
for systematic surveys for hostless SNe over the entire 
range of halo masses to study the spatial extent, stellar 
content and mass of diffuse light at all scales. 

We have recently completed the Multi-Epoch Nearby 
Cluster Survey (MENeaCS), one of whose many goals 
is to measure the SN la rate in the most massive, X- 
ray selected z ~ 0.1 galaxy clusters. In this paper, we 
report the discovery of four likely IC SN la which we 
have spectroscopically confirmed from among a sample 
of twenty-three in the MENeaCS survey. We use these 
to constrain the IC stellar population and stellar mass 
fraction in the MENeaCS cluster sample, along with the 
radial distribution of IC stars. The definition of ICL is 
necessarily ambiguo us, given that gala xies have no well- 
defined 'edge' (e.g. 
the point of view of 



Aba di et al.l 12006ft We thus adopt 



Rudick et al.l (|2010( ) and others: the 



ICL includes all low luminosity components of galaxy 
clusters, whether it is large scale diffuse luminosity, a 
tidal stream, or even the extreme outskirts of a cluster 
member. Where necessary, we adopt a flat cosmology 
with H =70 km s^Mpc^ 1 , O m =0.3, and fi A =0.7. 

2. SURVEY DESCRIPTION 

While a full description of the MENeaCS survey, along 
with our cluster SNe rates, will be presented in a sepa- 



rate paper (Sand et al. in prep), we briefly describe the 
relevant aspects of the survey here. The principal time- 
domain goal of MENeaCS is to measure the cluster SN la 
rate at z ~ 0.1, including both hosted and hostless, IC 
events. We take advantage of CFHT's queue schedul- 
ing to observe ~30 clusters per queue run (with 2 x 120 
s g' and r' images), and a total of 58 unique clusters 
throughout the year. These clusters were X-ray selected, 
and have 0.05 < z < 0.15. The monthly CFHT imaging 
is supplemented by follow-up imaging with the 2.3 m 
Bok telescope on Ki tt Peak with its 90Prime imager 
([Williams et al.ll2004l ). Viable cluster SNe la are identi- 
fied based on brightness and color selection, and spectro- 
scopically confirmed via monthly spectroscopy, primarily 
with the MMT. Target of Opportunity (ToO) observa- 
tions with Gemini were acquired specifically to follow up 
our hostless SN candidates, which were pursued without 
any prior cut on brightness or color. The survey began in 
February 2008 and spectroscopic observations ended in 
December 2009. The last of the imaging data was taken 
in 2010A in order to get the image depth necessary for 
the other scientific aspects of the MENeaCS program. 

During the course of the survey, we maintained deep- 
stacked images of all of our survey fields - consisting of 
all previous images at that position - for the purposes 
of faint host identification. For every SN candidate we 
attempt to identify a host galaxy using th e method of 
the SN Legacy Survey ([Sullivan et al.1 [20061 ). which uses 
a host-SN elliptical separation normalized by the candi- 
date host galaxy size to assign hosts - the dimensionless 
R parameter described in Eqn. 1 of Sullivan et al. (2006). 
Each host is characterized by an elliptical shape, as out- 
put by SExtractor, with semimajor and scmiminor axes 
(va and tb) along with a position angle, 6. We note here 
that Eqn. 1 of Sullivan et al. (2006) had a typographi- 
cal error. The corrected version of the formula - where 
(£sn,Z/sn) is the SN pixel position and (x ga \,y ga i) is the 
host center - should read 



R 2 



(i) 



2/SN - 2/gal, C xx — 

yy =s\n 2 {9)/r\ +cos 2 (0)/r!, 



where x r = x S n - x ga \, y r 
cos 2 (8)/r 2 A + sm 2 (9)/r 2 B , C, 
and C xy = 2cos(6») sin(0)(l/r| - l/r|). 

Unfortunately, the original typographic error propa- 
gated into our survey code, although as we shall see from 
Figure 1, it did not ultimately effect our ability to dis- 
tinguish between hosted and hostless SNe. Typically, 
the apparent detectable edge of a galaxy is at R ~ 3, al- 
though this is clearly a depth dependent statement. As in 
the SNLS, we initially classify a SN candidate as hostless 
if it has no potential hosts with R < 5. These candidate 
IC SNe are flagged and a visual search for a faint host 
is performed with our deep-stacked image, in case it was 
missed by our automated analysis. Very occasionally, the 
visual search revealed that a SN candidate with R > 5 
did seem to be associated with the outer halo of a galaxy 
(in retrospect, this was due to the typographical error 
mentioned above), and the SN was followed up during 
our classical spectropscopy observations. The individual 
discovery g' and r' band images are blinked and exam- 
ined by a human to verify that the IC SN candidate 
is not a slowly moving object. If the candidate is still 
viable, we triggered our Gemini ToO program, or oth- 



Intracluster Supernovae MENeaCS 



3 



erwise placed the hostless SN candidate into our 'high 
priority' spectroscopic follow up list. Given our deep 
imaging (see § I3.3|) . our adopted definition of a hostless 
SN should be comparable to or st r icter than recent IC SN 
work, such as iMcGee fe Baloghl (|2010( ) , who used twice 
the radius at which 90% of the light was enclosed (down 
to SDSS depths) as their working definition of a host- 
less SN. We also point out that, given the stacking of 
our images prior to SN discovery (with our shallowest 
stacks reaching M g ~ —12.5 mag; M r ~ —13.0 mag- 
see ?I4.2|). our host galax y limits are stricter than that of 

IMcGee fe B aloghl (l2010l) , which wer e at best M r 15.0 

mag, and similar to iGal-Yam et al.l (|2003l) . 

In general, our classically scheduled spectroscopic fol- 
low up was focused on likely cluster SN la. If a SN can- 
didate was newly discovered, its g — r color and mag- 
nitude had to be consistent with a SN la within ~30 
days of explosion at the cluster redshift, which ranged 
between -0.5 < g - r < 0.8 and 17.5 < g < 22.5 
(we have cali brat ed all survey imaging to the SDSS fil- 
ter system; § 13.11) and had to have a clustercentric dis- 
tance of R < 1.2i?2oo (we have derived -R200 for each of 
th e MENeaCS clusters u s ing th e X-ray scaling relations 
of iReiprich fe Bohringerl {2002), see 5 14.31) . Although 
every effort was made to spectroscopically confirm vi- 
able SNe la at higher clustercentric radius, this was not 
always possible. We also attempt to get spectroscopic 
confirmation of SNe even if they appear over 1 month 
old, based on our cadence, weather on previous spec- 
troscopic runs, and the light curve of the SN candidate. 
Note that all viable cluster SN la candidates, as defined 
above, were observed modulo telescope scheduling and 
weather. Other SN candidates (i.e. those that were not 
consistent with expected cluster SN la's colors or lumi- 
nosity, but still with g < 22.5) were followed up with 
lower priority, but this yielded 7 cluster core collapse 
SNe. These allow us to set a limit on the number of 
hostless to hosted core collapse SN e in the cluster en- 
vironment for the first time (§ 14. 3[) . a notion particu- 
larly relevant now that some st ar formation has been 
observed in the IC env i ronment dSun et aDl2007l 120101 : 
iSivanandam et al1l2010t iSmith et al.ll2010D~and has e ven 
been suggested by simulations (|Puchwein et al.ll2010D . 

MENeaCS has spectroscopically confirmed twenty- 
three cluster SN la (all within 3000 km s _1 of the cluster 
redshift), four of which have no host upon close inspec- 
tion. We plot the SN-to-host distance for our cluster SN 
la sample in Figure [IJ in terms of their effective radius, 
i? e , after host assignment via the elliptical radius pa- 
rameter as described above. The effective radius of each 
host galaxy is e stimated via the FLUX_R ADIUS parameter 
in SExtractor (|Bertin fc Arnoutsl 119961 ) (this parameter 
was set so that it reported the radius at which 50% of the 
galaxy light was enclosed, the traditional definition of ef- 
fective radius). The distances, in effective radii, of our 
four hostless cluster SNe la with respect to the nearest 
source in our SExtractor catalogs are shown as dashed 
lines (see § 13.31 for our limits on faint hosts associated 
with our IC SN candidates). Note that the SExtractor 
flux_radius parameter does not account for the point 
spread function (PSF) width, and so our R e values will 
be overestimates, with the corresponding host - SN dis- 
tance becoming larger. This plot reinforces the distinct 
nature of our IC SN candidates. 



3. OBSERVATIONS & ANALYSIS 

In this section, we will present the observations of our 
IC SNe along with our evidence for them not being re- 
lated to a visible host. In all, we have identified three 
cluster SN la which are excellent candidates for having 
progenitors which are IC stars, along with an intrigu- 
ing SN which may be associated with debris or other 
stripped material from a nearby cluster galaxy. We sum- 
marize the observations of each SN in Figures 2-5. The 
results of our spectroscopic analysis are summarized in 
Table [TJ and our magnitude limits for faint hosts are in 
Table H 

3.1. Discovery and follow up photometry 

We have implemented a nearly real time analysis sys- 
tem for discovering and flux calibrating SN candidates 
in our CFHT imaging. This pipeline, as discussed in 
§ also attempts to assign a host to every discovered 
transient. If no host is found, and this is confirmed via 
visual inspection of the raw and processed data, the SN 
is moved to the top of our spectroscopic priority list, 
and our Gemini ToO was triggered, if available. Since 
our imaging consisted of not only our CFHT/Megacam 
data, but also 2.3 m Bok and occasional Gemini ToO ac- 
quisition images, we chose to calibrate our images to the 
Sloan Digital Sky Survey g,r photometric system, which 
in general required both a zeropoint and color term. 

3.2. Spectroscopic confirmation 

We used our Gemini Observatory ToO program to con- 
firm three of our candidate IC event s using the Gemin i- 
North Multi Object Spectrograph (jHook et al.1 l2004f ). 
Long-slit spectra were taken with the R400 grating cen- 
tered at 680 nm, with the GG455 order blocking filter, 
covering ~465-890 nm. The l'.'O slit was used and was 
oriented so that a relatively bright star was also placed 
in the slit, in case the SN was not visible in acquisition 
images. In the case of AbelL399_3_14_0 (note that we 
adopt our internal MENeaCS identification for our SNe, 
although we list IAU circular identifiers when available 
in Table 1), the slit was oriented so that a redshift for 
the nearby cluster red sequence galaxy could be obtained 
simultaneously. All SN spectra were reduced in a stan- 
dard manner utilizing IRAFQ routines, including the 
Gemini-specific IRAF package. 

We have also u sed Hectospec on the MMT 
(|Fabricant et al.| [2005l in a queue mode to obtain a spec- 
trum of Abelll650_9_13_0. Hectospec is a 300 fiber (1'.'5 
diameter) spectrograph with a 1 square degree field of 
view. The 270 line grating was in place, giving wave- 
length coverage between ^360-800 nm. The spectrum 
was pipeline processed at the Harvard-Smithsonian Cen- 
ter for Astrophysics. 

The spectra of our IC SNe can be seen in Figures 2-5. 
We utiliz e the Supernova Identifi cation (SNID) software 
package (jBlondin fc Tonrvl l2"007') for classifying all of our 
spectra, which cross-correlates the input spectrum with 
a library of template SN spectra to both classify and de- 
termine the redshift of our SNe, based on the technique 

12 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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of lTonrv fc Davisl (fl979h. We mimic the four step proce- 
dure of Foley et al.l (|2009) for determining a SN's type, 
subtype, redshift and age using SNID. The reader is re- 
ferred to that work for the details of the process, and can 
see our results for the MENeaCS IC SNe in Table ED The 
best-matched SNID template spectrum is also listed in 
Table [TJ and is overplotted in Figures 2-5 for comparison 
with the observed spectrum. In general, the match be- 
tween template and spectrum is excellent, and we discuss 
the individual fits in § 13.41 

3.3. Deep stack images 

With the survey complete, we have stacked all CFHT 
imaging epochs without SN contamination to search for 
faint hosts in the vicinity of our IC SNe. In general, this 
consisted of imaging taken prior to SN discovery only; 
close examination of images taken a year or more af- 
ter SN discovery showed some residual SN light. Our 
final deep stack catalog for each field was made by me- 
dian combining these SN-free images, using the SWariEI 
software package, and then running SExtractor on the re- 
sulting image. Since the goal was to identify the faintest 
possible hosts, we set the SExtractor keywords De- 
tect_minarea to 5, Detect_thresh to 2, and Anal- 
YSIS_thresh to 2. No source was found within R < 5 
(see § [2]) confirming our original trigger on these events 
as potential IC SNe. We visually inspected the IC SN 
positions as well, again finding no host. 

To quantify the depth of these deep stack images, we 
implanted artificial stars (^50000 per image, over several 
iterations) with a PSF matched to that of the original 
image. We note that dwarf galaxies in the Coma cluster 
have a typical size of ~1 arcsecond (e.g. lKomiyama et al.l 
2002). Given that our cluster sample is at higher red- 
shift, and has a typical seeing of ~0. 8-1.0 arcseconds, 
our decision to quantify our sensitivity to faint hosts via 
a standard point source detection efficiency calculation 
is justified. The g magnitude of the artificial stars was 
drawn randomly from 21 to 28 mag, with uniform prob- 
ability. The g — r color was then randomly assigned over 
the range 0.5 < g — r < 1.5, again with uniform prob- 
ability. No artificial star was put within ten R e of any 
source in our original SExtractor catalog, as determined 
via the FLUX_RADIUS output parameter. Note that 
our statistics for sources >5i? e - more representative of 
our actual IC SN population as seen in Figure 1 - and 
>10i? e are identical (see § 4.1). After implantation, our 
detection efficiency as a function of magnitude was calcu- 
lated by re-running SExtractor on the image. We adopt 
the magnitude at which we recover 50% of our artificial 
sources as our image depth, and we report this in Ta- 
ble [5J after translating into absolute magnitudes given 
the redshift of each cluster. We inspected several images 
to confirm that the 50% recovery threshold in our SEx- 
tractor catalog was visually correct. We calculate what 
fraction of the clusters' luminosity is contained in dwarf 
galaxies below our detection limit, and the implications 
for that on the hostless nature of our IC SN candidates, 
in §0 

3.4. Notes on individual objects 

13 version 2.15.7; http://terapix.iap.fr/soft/swarp 



Here we present additional details for each of the IC 
SN candidates. 

3.4.1. AbeU1650-9-13-0 

Abelll650_9_13_0 was discovered 2009 December 14 
(UT dates are used throughout this paper) with g = 18.3 
mag and r — 18.6 mag, making it a prime target for 
spectroscopic follow up. A spectrum was taken with 
MMT/Hcctospec on 2009 December 21, which indicates 
a normal SN la near maximum light (Figure [2]). 

The SN position is -470 kpc East of Abell 1650's 
BCG. Deep stack images of Abell 1650 with no SN con- 
tamination indicate a depth of M g ~ —12.5 mag and 
M r ~ —13.0 mag at the cluster redshift (z — 0.0838). 
The SN is ~10 effective radii from the nearest object in 
our deep Abell 1650 catalog. 

3.4.2. Abell2495-5-13.0 

Abell2495_5_13_0 was discovered on 2009 May 23 with 
g — 23.0 mag and r = 22.4 mag, fainter than we typically 
followed up SNe in MENeaCS, but due to its lack of an 
apparent host and the fact that Abell 2495 had not been 
observed since October 2008, we triggered our Gemini 
ToO program. Further imaging on June 17 showed that 
the SN was in decline, just prior to our spectroscopic 
observations on June 18. The SN spectrum indicates 
that Abell2495_5_13_0 is a normal SN la roughly three 
months past maximum light (Figure [3]). 

The SN position is ~150 kpc from the BCG in Abell 
2495, and is the most centrally located in our sample 
(Figure [3]). The deep stack images indicate that there is 
no host galaxy at the location of Abell2495_5_13_0 to a 
magnitude of M g ~ —11.7 mag and M r ~ —12.4 mag 
at the distance to Abell 2495 (z = 0.0775). The nearest 
object in our Abell2495 deep catalog implies that the SN 
was over ~20 effective radii away. 

3.4.3. Abell85.6.08.0 

Abell85_6.08.0 was discovered on June 18, 2009 with 
g = 20.6 and r = 19.8. Subsequent photometry was 
obtained with both the CFHT and the Steward 2.3 m, 
indicating a SN well past maximum light. 

The Gemini/GMOS SN spectrum is shown in FigureSJ 
The best matched SNID template for Abell85_6_08.0 
is SN 1991bg, roughly five weeks past maximum, and 
the top six SNID matches to the spectrum are also 
of the SN 1991bg subtype. We are confident that 
Abell85_6_08_0 is an underluminous SN la based on OI 
and Ca II in absorption - which is typical for under- 
luminous SN la roughly 4-5 weeks past maximum (e.g. 
iTaubenberger et al.ll2008ft . We do not see the narrow Na 
I D line near 5700 A, typical of underluminous SN la. 
The SN type is supported by the declining light curve 
at the time the spectrum was taken, although no data is 
available near maximum light. If an underluminous SN 
la, Abell85_6_08 is likely ass ociated with an old stellar 
population (e.g. Howell 2001), which we briefly discuss 
in S OI 

Abell85_6.08.0 is ~600 kpc southeast of Abell 85's 
BCG. No host is detected down to M g ~ —11.2 mag 
and M r ~ —11.7 mag in our deep stack catalogs at the 
cluster redshift (z = 0.0578). The SN is -5 effective radii 
from the nearest object in our deep Abell 85 catalog. 
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3.4.4. Abell399-3M-0 

Abell399_3.14_0 was discovered November 6, 2008 
(preimaging with the 2.3 m Bok telscope on October 25, 
2008 also caught the SN on the rise). As can be seen 
in the left panel of Figure [5j the SN was relatively near 
a cluster red sequence member, which we will refer to 
as Galaxy 1. Nonetheless, our deep stack image at the 
time indicated that the SN was —7 R e from the poten- 
tial host's center, and so we triggered our Gemini ToO 
observations. Further analysis of the field indicates that 
the SN did occur > 7R e from the galaxy center based on 
SExtractor's PHOT_fluxfrac parameter. 

To investigate further, we used GALFIT (|Peng et al.l 
2002) to fit parameterized models to Galaxy 1. Several 
simple models with either a single de Vaucouleurs or Ser- 
sic profile, along with a central point source, were fit to 
Galaxy 1, all suggesting that the SN was > 5.5 effec- 
tive radii from the galaxy center. However, these models 
had significant residuals, suggesting that an additional 
component was necessary to fit the data. After some in- 
vestigation, the best model consists of two Sersic compo- 
nents, along with a point source at the galaxy's center. 
The larger Sersic component has an effective radius of 
2.2 arcseconds (3.3 kpc), suggesting that Abell399_3_14_0 
was only —3.8 effective radii from the galaxy center, at 
least with respect to this largest scale component. Given 
the Sersic index of this component, n = 7.7, ~17% of the 
light lies at larger radii than the SN radius. In Figure HJ 
middle panel, we show the two-component Sersic model 
of Galaxy 1 subtracted from the image. Note that there 
are low surface brightness features projected at large 
radii around Galaxy 1, possibly a plume or tidal debris, 
suggesting that it has recently interacted with another 
galaxy. Signs of disturbanc e and tidal features are com- 
mon in both field ellipticals (jvan Dokkumll2005HTal et al 



2009), and their cluster counterparts ([Janowiecki et al 



2010) 



To corroborate our GALFIT results, we implemented 
a non-parametric approach to measuring the light pro- 
file of Galaxy 1. The data around Galaxy 1 were binned 
into 1 pixel-wide annuli, with deviant pixels flagged and 
removed through an iterative procedure. A sky annulus 
was taken, 20-pixels wide, and this sky value was sub- 
tracted from Galaxy l's light profile. This allowed for a 
direct measurement of the fraction of light enclosed as 
a function of radius, with ~10% of Galaxy l's light at 
radii larger than that of the SN. This rough agreement 
between the GALFIT and non-parametric analyses indi- 
cate that Abell399_3_14_0 lies in the outskirts of Galaxy 
1 in projection. 

Given the association of Abell399_3_14_0 with Galaxy 
l's outer halo, and the appearance that Galaxy 1 has 
recently had a significant interaction, it is interesting to 
look at the relative velocities of Galaxy 1 (which was 
placed in the Gemini slit) and the SN. Precision red- 
shifts of SNe are difficult to measure due to their broad, 
evolving spectral features. That said, taking the me- 
dian redshift of the good correlations found by SNID - 
based on 20 normal SN la templates - gives a redshift of 
z = 0.0603 ±0.0012 (where the uncertainty is the stan- 
dard deviation of the good correlation redshifts) . 

SNID also has several galaxy templates that it can 
use to quantify galaxy redshifts, but since it smooths 



the data, we choose instead to use the IRAF task 
RVSAO (|Kurtz fc M ink! 1 19981) . which is bas ed on the 
cross-correlation algorithm of iTonrv fc Davisl (fl979l) . to 
precisely measure the redshift of Galaxy 1 (note that 
SNID gives consistent results with what follows, albeit 
with higher error bars). Several templates were used, 
based on absorption line spectra of nearby galaxies along 
with a template made with K-giant spectra in M31 glob- 
ular clusters. All of these templates give nearly identi- 
cal results, to within 5z = 0.0001. Takin g the template 
which gave the best r-statistic value (jTonry fe Davisl 
[1971 . we measure a redshift of z = 0.06439 ± 0.00010 
for Galaxy 1. 

While the velocity difference between Abell399_3_14_0 
and Galaxy 1 are formally more than 3-er apart, it is also 
insightful to plot the velocity difference between Galaxy 1 
and each of the good SN template matches which SNID 
found, which we do in Figure |HJ Note that all of the 
good SN template matches indicate a velocity difference 
>500 km s -1 . This velocity difference suggests that 
Abell399_3_14_0 is unbound to Galaxy 1. Assuming that 
Galaxy l's half light radius is between — 1 kpc (our GAL- 
FIT result based on a sing le i? 1 / 4 profile fit) and -2 kpc 
(based on our nonparametric analysis of the light profile), 
then we would expect Galaxy 1 to have a velocity dis - 
persion -100-200 km s" 1 (e.g. Ivan der Wellt~aT][200l . 
Using this, and the proje cted mass profile scaling rela- 
tion of lBolton et al.l (|2008f) based on strong gravitational 
lensing (which we deproject assuming spherical symme- 
try), we find that Galaxy 1 would have an escape velocity 
of -280-480 km/s at the radius of Abell399_3_14_0, pre- 
suming this is the true radius of the event. We grant 
that this is merely an e stimate and that th e application 
of the scaling relation of lBolton et al.l (|2008l) requires ex- 
trapolation, but it is indicative that Abell399_3_14_0 is 
at best loosely bound to Galaxy 1. 

We note that the redshift of Abell399.3_14_0 derived 
via SNID puts the SN just within our 3000 km s" 1 
velocity cutoff for designating cluster SNe. To inves- 
tigate if this velocity offset is reasonable, we searched 
the NASA Extragalactic Database (NED) for redshifts 
in the literature associated with Abell 399. A standard 
biweight estimator yields a cluster velocity dispersion of 
a = 1230^ 70 km s" 1 from 123 cluster redshifts (calcu- 
lated as in I Just et all 120101 ) . putting Abell399_3_14_0 in 
the 'outskirts' of the cluster in velocity space - but still 
within —2-3(7 given the SN's redshift uncertainty. With 
twenty-three cluster SN la in our sample, it is not unrea- 
sonable to have a SN with such a velocity offset from the 
cluster. 

The status of Abell399_3_14_0 leans toward it being an 
IC SN. In projection, the SN lies in the outer regions of 
Galaxy 1, with —10-20% of the galaxy's light at larger 
radii than the SN. On the other hand, Galaxy 1 appears 
to have associated tidal debris and the velocity offset 
between Galaxy 1 and the SN is larger than a simple 
estimate of the escape velocity at the projected radius. 
Indeed, finding a true IC SNe projected onto a galaxy is 
expected for samples roughly the size of MENeaCS (see 
§ 14. ip . The bulk of the evidence leads us to consider 
Abell399_3_14_0 an IC object based on these last facts. 
However, if one wants to adopt the stricter definition 
that IC stars are only those which are in the population 
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gravitationally dominated by the cluster rather than a 
local galaxy, than Abell399_3_14_0 could be a transitional 
object. If we consider the progenitor of Abell399_3_14_0 
to be part of a 'transition' stellar population becoming 
unbound (or remaining loosely bound) from its galaxy 
of birth, it d oes not appear to be alone in the literature. 
iSmithl (|1981l) reported the discovery of SN 19801, a SNIa, 
to be projected onto a stellar bridge between NGC4374 
and NGC4406 in the Virgo cluster, suggesting a second 
'transition' population object. As we calculate the IC 
stellar mass fraction in § 14.41 we will present our results 
both with and without Abell399_3 _14_0 as an IC SN. 

4. DISCUSSION & RESULTS 

In this section we will calculate the IC stellar mass 
fraction in our sample of clusters based on the relative 
numbers of IC and hosted SNe (§ I4.4|) and discuss the 
relative spatial distribution of these objects (§ I4.5[) . First, 
though, we quantify several observational biases which 
cause a preference for the discovery and follow up of IC 
SNe in MENeaCS (§l4"T|). Following this, we show that it 
is unlikely that our IC SNe population are coming from 
dwarf galaxies in the cluster, below our host detection 
limit, based on the depth of our deep stack imaging fields 
and standard cluster galaxy luminosity functions (§ 14. 2j) . 
We discuss the IC stellar population, with respect to the 
cluster galaxy population as a whole in § 14. 31 with special 
emphasis on what our SN sample contributes. 

4.1. Hostless supernova detection bias 

There are two potential biases in using the relative 
numbers of hostless and hosted cluster SNe to derive the 
IC stellar mass fraction in our sample, which we attempt 
to quantify here. These assessments will be used in § 14.41 
for properly calculating our IC stellar mass fraction based 
on the relative numbers of hostless and hosted cluster SN 
la, and we summarize them in Table |3] 

First, it is simply more difficult to discover SNe within 
hosted galaxies than it is in a seemingly blank portion 
of the sky. This would serve to boost the true num- 
ber of hosted SNe relative to unhosted SNe. To address 
this, we have carried out a series of detection efficiency 
simulations by adding fake SN la into our MENeaCS im- 
ages with realistic absolute magnitude and stretch distri- 
bution s using the template light curves of iNugent et al.l 
(2002). One set of these fake SN la was implanted with a 
distribution following the galaxy light in the image (out 
to a maximum effective radius of 5i? e ), while the other 
was purposely placed at R > 10R e away from any poten- 
tial galaxy host in order to mimic a hostless population 
(note that while this does not correspond to our declared 
limit of > 5i? e for our IC SNe, the statistics between im- 
planted objects at 5R e and > 10i? e are identical). The 
images, with their fake SN la implanted, were then run 
through our detection pipeline and the statistics of the 
recovered SNe were tallied. Details of this process will 
be presented when we calculate the cluster SN la rate 
(Sand et al. in preparation), but we use the relative de- 
tection efficiency between hosted and hostless SNe here 
to quantify how much easier it was to discover hostless 
cluster SN la in the MENeaCS survey. The hosted to 
hostless SN detection efficiency is 0.91 when integrated 
over all magnitudes and host-centric radii, and so a factor 
of 1.0/0.91 will be applied to any value of the number of 



hosted cluster SN la, Nh os tedi in relation to the number 
of IC SN la, N IC . 

Second, there is a spectroscopic bias which always 
biases the number of hosted SN la to lower values. 
Throughout the course of the MENeaCS campaign, we 
made every effort possible to follow up all transients that 
appeared to be hostless, as long as this was reasonable 
to do within our Gemini ToO program or our classi- 
cally scheduled spectroscopic time. Given our limited, 
classical spectroscopic resources (which were effected by 
scheduling and weather) we only followed up viable clus- 
ter SN la candidates with g < 22.5 during that time. 
We split our spectroscopic bias into two categories, and 
address each in turn. First, there is what we call 'spec- 
troscopic availability bias' - during the period in which 
our Gemini ToO program was active we were always able 
to spectroscopically confirm hostless SN candidates, even 
if we had no classical spectroscopic time due to schedul- 
ing, weather or instrument problems. Since we were able 
to get a spectrum of all of our viable cluster SN la can- 
didates as long as there were no losses due to scheduling, 
weather, or technical problems, this bias is easy to quan- 
tify. It should simply be the ratio of time when both 
classical spectroscopic time and Gemini ToO for host- 
less candidates was available over the total amount of 
time that spectroscopic follow up was available for host- 
less SN candidates - we calculate this to be 0.91. This 
means that another factor of 1.0/0.91 will be applied to 
any value of Nested when we calculate our IC stellar 
mass fraction. As a sanity check, another way of quan- 
tifying the spectroscopic availability bias is to directly 
investigate our transient database for good cluster SN la 
candidates (as described in § [2]) that were not followed 
up due to poor weather or scheduling; three such objects 
were found. Based on the 0.91 fraction of classical spec- 
troscopic coverage found above, we would expect one or 
two such objects - perfectly consistent, especially given 
that good SN candidates occasionally turned out to be 
background/foreground SNe or even core collapse SNe. 

The second type of spectroscopic bias is a 'hostless 
follow up' bias - we were willing to spectroscopically 
confirm fainter, apparently hostless, SN candidates with 
Gemini but would avoid pursuing similar hosted candi- 
dates with our classical spectroscopic time. The simplest 
way of addressing this hostless follow up bias is to strike 
Abell2495_5_13_0 from our calculations of the IC stel- 
lar mass fraction (§ I4.4|) . With a discovery magnitude of 
g = 23.0, r = 22.4, and declining light curve, such a clus- 
ter SN la caught so long after maximum would not have 
been spectroscopically followed up and confirmed if not 
for our Gemini ToO time. All other IC SNe would have 
been high priority targets in our classical spectroscopy 
queue due to their magnitude and color. 

There is another bias, not considered in previous work, 
due to true IC SNe being projected onto a false host 
galaxy, whether it be a foreground/background galaxy 
or a cluster member. To gauge the potential magnitude 
of this 'filling factor' effect we calculated the fractional 
image area enclosed within the dimensionless elliptical 
radius R < 5 (see § [2]) of any potential host in a repre- 
sentative sample of clusters, finding that on average ~10 
to 15% of the image area is covered by a potential 'host 
galaxy'. Additionally, when going from the cluster out- 
skirts to the inner regions of the cluster we find that the 
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filled area increases by ~ 2 — 3%. Of course, this fill- 
ing factor estimation method assumes that IC SNe are 
randomly distributed in the cluster, when in reality they 
could be more strongly associated with galaxies if they 
have been recently stripped; we do not quantify this plau- 
sible scenario further. Our best estimate suggests that 
<1 of our apparently hosted SN may be an actual IC 
SN projected onto a galaxy's position. One way to miti- 
gate this effect is to carefully compare SN and potential 
host redshifts (although this is limited due to the broad, 
evolving features in SN spectra), which we have done 
for the entire hosted SN set in the MENeaCS sample, 
of which ~30% have a measured host redshift (Sand et 
al. in prep). Besides Abell399_3_14_0 and its large ve- 
locity offset with its possible host galaxy (§ 13.4.41) . which 
we consider the first object in this class, we have found 
no SN with significant velocity offsets from their puta- 
tive host. Given the case of Abell399_3_14_0, we apply 
no additional correction factor since the expected size of 
the effect is comparable in size to that observed. We 
gauge its impact by presenting our results both with and 
without it in the IC sample. Future work, in the era of 
large time domain surveys, will likely turn up more such 
objects. 

To summarize, we will apply three correction factors, 
based on observational biases, to our IC stellar mass frac- 
tion analysis in § 14.41 which we show in Table [3] First, 
we strike Abell2495_5_13_0 from our IC stellar mass frac- 
tion calculations for being too faint in comparison to our 
hosted SN la. We then apply two separate factors of 
1.0/0.91 whenever a Nested term is necessary: once for 
the relative difficulty of detecting hosted SNe in compar- 
ison to hostless events, and once for the relative spectro- 
scopic availability of our IC SNe target of opportunity 
observations and our classically scheduled spectroscopic 
time for which we followed up our hosted SNe. 

4.2. Host limits and the cluster luminosity function 

Given our host detection limits and a cluster galaxy 
luminosity function, one can determine the fraction of the 
cluster's luminosity contained in dwarfs fainter than our 
observational limits. From this, the probability of having 
the number of apparently hostless SNe in the MENeaCS 
survey can be assessed, assuming that the SN la rate is 
roughly proportional to luminosity. 

As a starting point, we choose the luminosity func- 
tion and normaliza t ion of the Virgo Cluster measured by 
iTrentham fc Tullvl (|2002D . which has a faint-end slope 
of a — —1.03, and was measured down to Mr = —10 
mag. By integrating their luminosity function up to 
our detection limit for faint hosts, and normalizing it 
by the total luminosity, the fraction of cluster luminos- 
ity missed is calculated (not including the contribution 
from IC stars), and we report this value for each IC field 
in Tabic Note that in converting our reported M r 
magnit udes in Table [21 to the M r luminosity function 
used bv lTrentham fc Tulrv|f 2002), we made use of the fil- 
ter conversion relations of iBlanton fc Roweisl (j2007l ) and 
the typical color s of dw arf galaxies in SDSS reported by 
lEisenstein et all (|2001ft . 

A variety of cluster luminosity function faint end slopes 
have been measured in the literatur e - for a recent com - 
pilation for the Coma Cluster, see Mil ne et al.l ([2007) . 
The results of the last ten years indicate that the Coma 



Cluster has a faint end slope between —1.5 < a < —1.0. 
We use a value of a = —1.5 to give a realistic upper limit 
on the fraction of the cluster luminosity in faint galaxies 
below our deep stack d etection limits, aga i n usi ng the 
normalization terms of ITrentham fc Tullvl ([2002). We 
list these values in Table [2j 

Assuming that the SN rate is proportional to luminos- 
ity, it is unlikely that any of our IC SN candidates are ac- 
tually hosted by faint dwarfs below our detection limits. 
If the typical faint end slope were a = —1.0 3, as mea- 
sured for Virgo bv ITrentham fc Tullvl (|2002l ). the typi- 
cal amount of galaxy light below our detection threshold 
is 0.04%-0.12% of the total, suggesting that we would 
expect < 0.03 SNe below our detection limit. Even 
given our upper limit scenario, with a faint end slope of 
a = —1.5 and our shallowest detection limit implies that 
we would expect < 0.4 SNe la hosted by galaxies below 
our detection threshold. A faint end slope of a ~ —1.85 
is necessary for it to be plausible that all of our IC SN 
candidates were hosted by dwarf galaxies below our de- 
tection limit. 

We note that at some level all of these extrapolations 
are unphysical, and should only be taken as indicative. 
As is well known, a Schechter function with a faint end 
slope steeper than a = —1.0 has a diverging number 
of dwarf galaxies, while faint end slopes steeper than 
a = —2.0 diverge in luminosity. The faint end lumi- 
nosity function in clusters must truncate at the faintest 
luminosities, strengthening our conclusion that our des- 
ignated SNe are truly IC events. 

Of course, these strong constraints on faint dwarf 
galaxy hosts implicitly assume that the SN la rate is pro- 
portional to stellar luminosity at the faintest luminosi- 
ties. Recent results from the Lick Observatory Supernova 
Search (LOSS), however, suggest that there is a strong 
dependence on the SN la rate and the 'size' (and by ex- 
tension, luminosity) of a galaxy (|Li et al.l 120101 ) . This 
rate-size relation indicates that smaller galaxies have a 
higher SN rate per unit luminosity than larger galaxies. 
The physical mechanism for this effect is unknown, al- 
though it may be due to a h igher specific s tar formation 
rate among smaller galaxies ()Li et al.ll2010l ). If true, then 
the rate-size relation is not likely to hold in the cluster 
environment, since dwarf galaxies in cluster cores have 
had their st ar formation shut off by environmental pro- 
cesses (e.g. iPennv et al.ll2010| ). We also note that the 
targeted search strategy of LOSS did not include galax- 
ies as faint as M r > —13 mag, as we are discussing 
here. Additionally, due to their extreme stellar den- 
sities and low metallicities, globular clusters may have 
an enhanced SN la rate per unit mass, perhaps as high 
as a factor of ~10 flShara fc Hur lev! 12002]: | Ivanova et al.1 
[200a[Ross"wog et al.ll2008t iPfahl et al.ll2009f ) - intraclus- 
ter globular clusters have be en observed in the Virgo 
Cluster (e.g. iLee et all 120101) . Any effect of the rate- 
size relation on dwarf galaxies below our detection limit 
would push our true IC stellar mass fraction to lower 
values, as would any enhancement of the SN la rate in 
globular clusters. These possibilities make the case for 
even deeper high-resolution imaging in the vicinity of our 
IC SN candidates very exciting. 
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4.3. What is the make-up of the intracluster stellar 
population? 

A basic question when trying to use IC SNe (or for that 
matter, any tracer of IC stars) as a measure of the IC 
stellar mass fraction is the respective stellar populations 
of the cluster galaxies and the IC stars. Much recent 
work has pointed out that the SN la rate of a stellar 
population may have two components; one proportional 
to the ste llar mass, and the other to the star formation 
rate (e.g. iMannucci et al J 120051 : iScannapieco fc Bildstenl 
120051) . Alternatively, there may be a continuum of de- 
lay times between the birth of a st ellar population and 
when SN la expl osions occur (e.g. iPritchet et all l200l 
iMaoz etaD 120101 ). If, for instance, the IC stellar popu- 
lation was composed of solely old stars, while the cluster 
galaxy population had a mix, including continuing small 
levels of star formation, then the simple fraction of IC 
SNe out of the total would not be a good measure of 
the IC stellar mass fraction. Both the observational and 
numerical picture of the IC stellar population are varied 
and difficult to interpret, and we attempt to summarize 
the current state of knowledge here. 

Observationally, the most direct view of the ICL 
population has been in Virgo via Hubble Space Tele- 
scope imaging of i ts red giant branch population 
(|Williams et al.1 120071 ). Using color magnitu de diagram 
matching techniques, iWilliams et al.l (|2007l ) concluded 
that the ICL population in Virgo (at least in the field 
studied) is mostly old (>10 Gyr) with a wide range 
in metallicity (-2.3 < [Fe/H] < 0.0) but also has a 
young, relatively more metal-rich component (<10 Gyr; 
[Fe/H]> 0.5) as well. As a caveat, Virgo has a mass of 
-1.7xlO 14 M dRines fc Diaferiol 120061) . a nd is dynam- 
ically unrelaxed (e.g. iBinggeli et al.l 119851) . and so may 
not have ICL properties analogous to the MENeaCS 
cluster sample, which has a median cluster mass ^3.8 
times higher (using t he X-ray scaling relation fo r the 
flux-limited sample of Rciprich & Bohringcr 2002, and 
their bisector regression results). There is also observa- 
tional evidence for current star formation in IC space; 
the cold interstellar medium, removed from its galaxy 
due to ram - pressure stripping, can form stars (|Sun et all 
I2007L 120101 ISivanandam et al.llMot iSmith et afl T2010). 
Broad band color imaging of the ICL, on the whole, 
indicates that the ICL and the cluster galaxies have 
similar stellar populations, with some studies finding 
colors that are more similar to cluster galaxy out- 
skirts, while others finding colors as red as the interi- 
ors of early type galaxies; given the large uncertainties, 
these color measurements are br oadly consistent with 
the cluster galaxy population (e.g. iGonzalez et al1l2000l: 
Zibetti et all 120051: iKrick et all 120061: IKrick fc Bernstein! 
20071: iDa Rocha et al.l 120081: iRudick et al.l 120101 among 
others). 

Simulations have shown a similar range in results. For 
instan ce, si mulations of IC L build up by Wi llman et al.l 
(|2004l) and iMurante et~al"1 (|2004l) have suggested that 
IC stars are more metal poor and older than the gen- 
eral cluster galaxy population on average, respectfully. 
Other simulations indicate that IC stars have very sim- 
ilar properties to those in clu ster early type galax- 
ies (|Sommer-Larsen et al.l [2005) . Intri guingly, and in 
line with some recent observations (|Sun "etaLl [20101: 



ISmith et al.l 120101 ) . iPuchwein et~aTI (|2010l ) suggest that 
there may be some 'intracluster star formation' coming 
from gas stripped from small infalling halos. All of the 
simulation results suffer from so me overcooling, includ - 
ing those with AGN feedback (e.g. IPuchwein et al.ll2010f ). 
and can not reproduce the observed fraction of stars in 
cluster galaxies. Thus simulation results should be taken 
as only indicative. 

A crude hint of the make-up of the ICL population can 
be gleaned from the MENeaCS SN sample. One of our IC 
SNe, Abell85_6_08_0, has been identified as a SN1991bg- 
like event, while none of our hosted SN la can be placed in 
this class (although our sparse light curve sampling and 
single spectrum per SN may preclude precise classifica- 
tion). iLi et all (|2010l) report that ~18% of all SNe la are 
of the SN1991bg-like variety in a volume-limited survey, 
most appropriat e for our clust er SN search. They prefer 
early type hosts (|Howellll20011 ) , and thus presumably are 
associated with an old stellar population (although they 
may have a distribution of delay times). Note that we 
can not distinguish between scenarios in which the stel- 
lar population was stripped while still young - evolving 
passively in IC space since - or if the stellar population 
was stripped more recently with the old stars developing 
in a galaxy evironment. Given the small SN numbers 
involved, the most robust conclusion is that there is an 
old stellar population component to the IC light. Fu- 
ture large time domain surveys may be able to indirectly 
probe the stellar population of the ICL via the mix of 
their SNe population. 

As an aside, it should be mentioned that no IC core 
collapse SNe were discovered in MENeaCS, while seven 
such hosted events were. While the focus of MENeaCS 
was on SNe la, we agressively followed up all of our IC 
SN candidates. In this situation, where hosted candidate 
core collapse events were given lower spectroscopic pri- 
ority but not their hostless counterparts, we can place 
an upper limit of fic.cc < 0.35 (95% confidence upper 
limit) on the fraction of IC core collapse S Ne. This num - 
ber was derived using binomial statistics (Gehrels 1986), 
whose applicability is briefly discussed in § 14.41 If the 
core collapse SN rate is proportional to the current star 
formation rate (which it should be if core collapse SN 
progenitor s are from massive, short lived stars - for a 
review see IGal-Yam et al.l I2007D , then this fraction rep- 
resents the relative star formation rate in the IC popu- 
lation to that hosted by the cluster galaxies. While this 
constraint is not stringent, we again emphasize that fu- 
ture work associated with the large time domain surveys 
will be able to use IC SNe to probe the stellar population 
of the ICL. 

4.4. The intracluster stellar mass fraction 

The simplest assumption we can make is that the IC 
stellar population is similar to that of the cluster galax- 
ies. Our presentation of the IC stellar mass fraction in 
our sample will then focus on the fraction of IC SNe to 
the total (corrected, of course, for the biases discussed 
in § 14. 1|) as the proper tracer. Nonetheless, we will also 
present a conservative upper limit in case the IC stellar 
population is more s imilar to that of the older, early type 
clust er galaxies (e.g. lWillman et aLll2004J : IMurante et al.l 
120041 ). by utilizing those cluster SN la hosted by red se- 
quence galaxies (see below) . As we describe below, bino- 
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mial statistics are used to derive confidence limits, given 
that intracluster and hosted SNe are each distinguish- 
able, random events. We implicitly assume that the SN 
la rate is proportional to luminosity, and any deviation 
from this among globular clusters and dwarf galaxies be- 
low our faint host detection limit will push our intraclus- 
ter stellar mass fraction to lower values, as discussed in 

In § 14.11 we quantified our known sources of bias associ- 
ated with discovery and confirmation of IC versus hosted 
cluster SN la (summarized in Table [3]). Accordingly, two 
factors of 1/0.91 must be applied whenever the number 
of hosted cluster SN la is compared with the number of 
IC SNe, and Abell2495_5_13_0 should be struck from the 
sample since it would not have been followed up dur- 
ing our normal, classically scheduled spectroscopic time. 
This leaves three IC SNe out of a total sample of twenty- 
two, keeping in mind that we will present results with 
Abell399_3_14_0 both in and out of the IC sample. If we 
wish to include only those cluster SNe within -R200 (which 
we prefer, given that we begin becoming spectroscopi- 
cally incomplete at R > 1.2i?2oo)j where we have derived 
i?200 for each of the MEN eaCS clusters again us i ng th e 
X-ray scaling relations of iReiprich fc Bohringerl {2002), 
then the total sample shrinks to sixteen. We note that 
all but six of the MENeaCS clusters have been searched 
for SNe out to i?2oo, given the CFHT/Megacam field of 
view. Even for these clusters, we were able to cover out 
to ~ 0.9i?2oo, and we thus apply no additional correction 
factors to our SN numbers. 

If the IC stellar population is exclusively old, we need 
a comparable SN number among the hosted population 
in order to constrain the IC stellar mass fraction. To 
assess how many of our hosted cluster SN la were as- 
sociated with old stellar regions, we have constructed 
red sequence catalogs for each of our clusters (Bildfell et 
al. in prep), with a typical scatter in the color magni- 
tude relation of ^0.03 magnitudes. Seven (eight if one 
includes Abell399_3_14_0 among the hosted SNe) of the 
cluster SN la hosts are consistent with the red sequence 
of their cluster (five are within i?20o! six if y° u include 
Abell399_3_14_0). We will take our red sequence galax- 
ies to represent an exclusively old stellar population, but 
acknowledge that cluster red sequence galaxies may have 
an elevated SN la r ate compared to those in the field 
(iSharon et all 120071 : IMannucci et all 120081: IPiTdav et al.l 
l2010a|) . perhaps due to residual star formation or a sub- 
tle metallicity effect. We will address this issue with 
MENeaCS data in future work, although all of the ap- 
propriate studies su ffer from small number statistics. As 
iRudnick et al.l (|2009l ) measured ~75% of the clusters' to- 
tal light originating in red sequence members in their 
sample of SDSS clusters at (z) — 0.07, and younger 
galaxies will have a higher mass to light ratio, we conser- 
vatively assume that ^75% of our clusters' stellar mass 
is in red sequence members in order to make a correc- 
tion to represent the clusters' old stellar population as a 
whole. We apply this additional correction factor, which 
has the form 1 /0.75 whenever the number of SN la hosted 
by cluster red sequence galaxies is used to calculate our 
limits on the IC stellar mass fraction. This correction fac- 
tor likely biases the IC stellar mass fraction high when 
calculating a conservative upper limit for that quantity. 

We list our IC stellar mass fraction, given our various 



scenarios, in Table 0] To be explicit, the IC stellar mass 
fraction, fic,*, is found using 



Nic + {Nhostedbobsfpop) 

where Nig is the number of IC SN la and Nested is the 
number of hosted SN la. The quantity b bs is the to- 
tal observational bias factor determined in § 14.11 and is 
1/(0. 91) 2 throughout this work. Finally, f pop is an ad- 
ditional correction factor representing the fraction of the 
parent stellar population that the hosted SN sample is 
probing. The value of f pop is unity when we consider the 
IC stellar population and that of the cluster galaxies to 
be identical, and f pop = 1/0.75 when we are calculating 
the IC stellar mass fraction assuming that the IC stellar 
population is exclusively old, and use the number of SN 
la hosted by red sequence galaxies to represent the hosted 
old stellar population in the clusters - see the previous 
paragraph for details. All of our confidence limits were 
derived based on binomial statistics, since we are measur- 
ing two different kinds of distinguishable, random events 
(i.e. intracluster a nd hoste d SNe) . In particular, we use 
the expressions of iGehrelsl (|1986l ) for finding confidence 
limits on the IC stellar mass fraction, fic,*- We present 
all of our IC stellar mass fraction numbers as 68% con- 
fidence limits, while the red sequence derived IC stellar 
mass fractions are 84% one-sided confidence limits. 

Focusing on our results within -R20O7 we find an ICL 
stellar mass fraction of O.lGlQ' Jg (O.IIIq;^ if you do not 
include Abcll399_3_14_0 as an IC SNe) . If we use our con- 
servative measure of the ICL stellar mass fraction (count- 
ing only our SNe la hosted by red sequence galaxies and 
including Abell399_3_14_0 among our IC SNe), we derive 
an upper limit on the IC stellar mass fraction of < 0.47, 
which we consider the absolute upper range that our re- 
sults allow. 

This IC stellar mass fraction is consistent with 
the bulk of the results in the literature, which 
generally find ~10-30% at the cluste r scale (e.g. 
I Gal-Yam et aLlUM IGonzalez et all [2001 IZibetti et all 
120051: IKrick fe Bernstein! 120071 ). IGonzalez et all (|2007f ) 
measured an ICL fraction versus cluster velocity disper- 
sion relation for a sample of 23 groups and clusters. We 
apply this scaling relation to our own sample to deter- 
mine our expected ICL fraction, correcting the Gonzalez 
et al. numbers for the available field of view (using Fig- 
ure 5 in Gonzalez et al. 2007), and the velocity dispersion 
of the clusters in our sample (co nverting from th e X-ray 
luminosity using the relation of lWu et all (fl999)). Ad- 
ditionally, we apply a rough correction fact o r of .8 to 
account for the fact that the IGonzalez et al.l ([2007D ICL 
scaling relation actually includes the BCG luminosity, 
rou ghly ~20% of the tot al BCG+ICL luminosity budget 
(see lGonzalez et aLll2005l Figure 7) . We find an expected 
ICL fraction of ^0.19, which is consistent with our cur- 
rent, SN-based measurement. 

This can be seen more clearly in Figure [71 where we 
compa re our SN IC re s ults w ith those at the group 
scale (jMcGee fc Balogh] 1 2 01 0T) and the d i rect s urface 
brightness measurements of Gonzale z et all ((2007). Two 
simple correction fac tors were made to the original 
IGonzalez et all (|2007l) numbers for the purposes of this 
plot. First, a factor of 0.8 was applied to account for the 
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fact that the Gonzalez et al. numbers were presented as 
a BCG + ICL fraction, rather than the ICL alone; the 
factor of 0.8 is the rough fraction that Gonzalez et al. 
found themselves as the typical fraction of the BCG + 
ICL belonging to the ICL alone. Second, we renormalize 
their ICL numbers so that they reflect the fraction within 
R < R200, rather than the listed R < -R500 value (using 
Figure 5 of Gonzalez et al. 2007). We present two num- 
bers for the McGee et al. intragroup SNe results. The 
lower, dashed box represents the raw McGee et al. mea- 
surement of 19 intragroup SNe out of a total sample of 
59. The upper, dotted box is the final intragroup stellar 
mass fraction they reported after assuming that the IC 
stars are universally old, the hosted group stars contained 
a mix of stellar ages, and then applying a correction fac- 
tor to their raw SN numbe rs using the t wo-com ponent 
SN la rate measurement of [Dildav et al.l (j2010bj ). Note 
that this should be analogous to our own estimation of 
the upper limit of the IC stellar mass fraction, assum- 
ing that the IC stars are universally old. The velocity 
dispersion range of the McGee et al. group sample was 
derived by taking their quoted mean group mass, with its 
sta ndard deviat i on, an d converting t o a using Equation 
6 of lYang et all (|2007[ ). Note that the I Yang etaLI (|2007t ) 
group catalog was the origin of the McGee et al. group 
sample. 

The MENeaCS box in Figure [7] represents our ICL 
stellar mass fraction result (within -R200) including 
Abell399_3_14_0 as an IC SN. The upper limit arrow rep- 
resents the limit on the IC stellar mass fraction assuming 
that all IC stars are from an old stellar population (while 
including Abell399_3_14_0 as an IC SNe). The combined 
SN results at the group scale with our own confirm the 
declining IC s tellar mass fract i on, as a function of cluster 
mass, seen by I Gonzalez et all (|2007l) . 

The decline in IC stellar mass fraction as a function of 
halo mass may point to the assembly history of massive 
clusters. Speaking in broad terms, if one assumes that 
groups are the dominant building block of clusters then 
the groups that made present day clusters could not have 
had the high IC stellar mass fraction that we see today, 
unless some mechanism could remove intragroup stars 
during the cluster formation process. Moderate redshift 
measurements of the intragroup stellar fraction should 
be lower than those seen at low redshift to fit in with 
expectations from hierarchical structure formation. 

4.5. The spatial distribution of intracluster supernovae 

In Figure [5] we present the distribution of SNe la as 
a function of projected distance from the cluster center, 
scaled by i?2oo- Represented by the black histogram, it 
is apparent that the IC SNe are more centrally concen- 
trated in the cluster than the general, hosted popula- 
tion (keep in mind the MENeaCS SN coverage begins 
to be incomplete at R > R200) ~ although this is sub- 
ject to small number statistics. The fact that the ICL is 
more centrally concentrated than the cluster's galax ies 
has been observed previously (e.g. iZibetti et al.ll2005f l. 

The central concentration of the IC SNe is not a clear 
cut marker for the underlying stellar population of the 
ICL. On the one hand, this may favor a scenario in which 
the IC stellar population is primarily old, since low red- 
shift cluster cores are dominated b y red sequence galaxie s 
with old stellar populations (e.g. lUrquhart et aLll2010l ). 



On the other hand, iSmith et al.l (|2010f ) have shown, in 
their study of ultraviolet tails in the Coma cluster, that 
it is those few blue galaxies in Coma's central regions 
(modulo projection effects) which are predominantly be- 
ing stripped and contributing young stars to IC space. 
Our IC SNe distribution can not distinguish between 
these two scenarios. 

In a precursor cluster SN survey to the current work, 
iSand et all (|2008t) found 3 IC SN candidates at R > i? 2 oo 
within their transient sample, and none at smaller radii, 
suggesting a relative deficit if the candidates proved au- 
thentic. Although it is plausible that infalling groups 
have their own detectable IC stellar population, this 
seems to be at odds with the current work, which has 
a similar radial coverage. Table[T]shows that our IC SNe 
were between ~150 and 600 kpc from the cluster cen- 
ter (as measured from the BCG position). We suggest 
that the IC SN ca ndidates at high clustercentric radii in 
ISand et all ([2008) were likely background events in need 
of spectroscopic confirmation, although we can not rule 
out that individual events were indeed associated with 
an extended IC stellar population. Indeed, although not 
discussed earlier, MENeaCS had one viable IC SN la 
candidate which was followed up with our Gemini ToO 
program which turned out to be a background SN la at 
z = 0.2. Deep imaging may yet uncover a faint host at its 
position, and other appare ntly hostless SNe are routinely 
reported in the literature ([Sullivan et al.|[2006|) . 

5. SUMMARY AND CONCLUSIONS 

The MENeaCS cluster SNe survey has spectroscopi- 
cally confirmed a sample of twenty-three cluster SN la, 
four of which have no hosts upon close inspection. We 
use this sample of hostless SNe to constrain the IC stel- 
lar mass fraction, and the extent of the ICL in clusters. 
We have presented confirmational spectroscopy of these 
IC SN candidates, concluding that Abell85_6_08_0 was 
likely a SN 1991bg-like event, suggesting that a compo- 
nent oftheIC_stel]ar_p_opulation is old, as seen by others 
(e.g. lWilliams et al.ll2007T ). Deep CFHT/Me gacam imag- 
ing indicated that one of our events, Abell399_3_14_0, is 
spatially coincident with the outskirts of a nearby red 
sequence galaxy, but the velocity offset between the SN 
and its putative host (along with evident tidal debris in 
the region) show that it is likely unbound and thus an IC 
SN. We characterized our detection limit for faint host 
galaxies by implanting artificial stars into our deep stack 
images and found detection limits of M g > —12.5 and 
M r > —13.0 in a ll cases. Using the Virgo Cluster lumi- 
nosity function of lTrentham k, Tuflvl (|2002D . we calculate 
that < 0.1% of the clusters' luminosity is in dwarf galax- 
ies below our detection limit. Assuming a very steep faint 
end slope of a = —1.5 increases this number to < 1.7%, 
making it very unlikely that more than one of our IC SN 
la candidates was hosted by a faint dwarf, as long as there 
is no mechanis m for boosting the SN la rate am ong very 
faint hosts fe.g lPfahl et~alll2009l: iLi et al.ll2010D . All four 
IC SNe occurred within ~600 kpc of the cluster center, 
while our hosted SNe are seen out to -R200 and beyond, 
illustrating that the ICL is more centrally concentrated 
than the cluster galaxies. 

After accounting for observational biases which made 
our IC SNe easier to discover and follow up, we calcu- 
late an IC stellar mass fraction of 0.16to'og (68% confi- 
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dence limit) for all objects within i?200; assuming that 
the cluster galaxies and the IC stars have a similar mix 
of stellar populations. If we instead assume that the IC 
stellar population is exclusively old, and that the clus- 
ter galaxies have a mix of stellar ages, we derive an up- 
per limit on the IC stellar mass fraction of < 0.47 (84% 
one-sided confidence limit). This is in excellent agree- 
ment with the majority of IC measurements in the lit- 
erature. These results, along with the IC stellar mass 
fr action at the group scale provided by the SNe study 
of IMcGee fc Baloghj (|2010t ). allow us to confirm the de- 
clining stell ar mass fraction as a f unction of halo mass 
reported by iGonzalez et al.l ()2007l ) . IC and intragroup 
studies at higher redshift will further constrain the as- 
sembly history of these massive halos. Indeed, one ex- 
citing avenue is the pursuit of IC SNe to higher red- 
shift, where cosmological surface brightness dimming will 
greatly reduce the utility of traditional direct measure- 
ments of the ICL. 

Future prospects for hostless SN studies are excellent. 
With the advent of large, dedicated time domain surveys 
such as the Palomar Transient Factory, SkyMapper, Pan- 
Starrs, La Silla-QUEST Variability Survey and the Large 
Synoptic Survey Telescope, the number of hostless SNe 
should increase rapidly. These additional numbers will 
provide the statistical confidence necessary to truly con- 
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Fig. 1. — A histogram of the SN to host galaxy distance for the cluster SNe la in the MENeaCS sample, expressed in units of the host's 
effective radius, R e . The four IC SNe discussed in the current work, shown as dashed lines, lie at a host - SN distance > 5i? e from the 
nearest possible host. 
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Fig. 2.— Left: Discovery r' image of Abelll650.9.13.0, which is located ~470 kpc East of the BCG in Abell 1650 (z = 0.0838). North 
is up and East is to the left. Middle: A zoom in of the SN location on our deep stack, SN-free image. No host is apparent down to 
M g = -12.47 mag, M r = -13.04 mag at the redshift of the cluster (Tabled. Right: MMT/Hectospec spectrum of Abelll650_9_13.0, 
along with the best- fitting SNID template, SN2002er, which is overplotted in red. The SN is a normal SN la at z = 0.0836 (see Table [TJ- 
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Fig. 3.— Left: Discovery r' image of Abell2495.5.13.0, which is located ~150 kpc Northeast of the BCG in Abell 2495 (z = 0.0775). 
Middle: A zoomed in view of the SN location in our deep stack, SN-free image. No host is apparent down to M g = —11.72, M r = —12.37 
mag. Right: Gemini/GMOS spectrum of Abell2495_5_13_0, along with the best-fitting SNID template, SN1998aq, which is overplotted in 
red. 
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Fig. 4.— Left: Discovery r' image of Abell85.6.08.0, which is located ~600 kpc Southeast of the BCG in Abell 85 (z = 0.0578). Middle: 
A close up view of the SN location in our deep stack, SN-free image. No host is apparent down to M g = — 11.15, M r = —11.68 mag. 
Right: Gemini/GMOS spectrum of Abell85_6_08_0, along with the best-fitting SNID template, S N91bg , which is overplotted in red. Even 
though the best-matched SNID spectrum is not as well fit as our other spectra, we argue in § 13.4.31 that Abell85_6_08_0 is likely of the 
underluminous, SN 1991bg subtype. 
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Fig. 5.— Left: Discovery r' image of Abell399.3.14_0, which is located ~620 kpc Southwest of the BCG in Abell 399 (z = 0.0718). Note 
the relatively nearby cluster galaxy, dubbed Galaxy 1. Middle: A close up view of the SN location in our deep stack, SN-free image, where 
the best GALFIT model of Galaxy 1 has been subtracted (along with several other nearby galaxies). Not ice the debris uncovered after 
the GALFIT model s ubtraction, reminiscent of that seen in local cluster early types ( Janowiccki ct al. 2010). We discuss the IC status of 
this object in 5 13.4.41 Right: Gemini/GMOS spectrum of Abell399.3.14_0, along with the best- fitting SNID template, SN2003du, a normal 
type-la, which is overplotted in red. 
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Fig. 6. — A histogram of the velocity offset between Abell399_3_14_0 and Galaxy 1 for the twenty good SNID template matches to the 
SN. None of the template redshifts are within 500 km s _1 of the galaxy. We argue that Abell399_3_14_0 is a part of a stellar population 
which either is unbound or is becoming IC stars since it lies in the outskirts of Galaxy 1, has a significant velocity differ ence w ith that 
system, and Galaxy 1 itself has associated low surface brightness plumes indicating a recent interaction (see discussion in § I3.4.4[l . 
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Fig. 7. — A comparison of the IC stellar mass fraction measured in the current work (solid lined box for our numbe rs at R < -R200 
includ ing Abell399_3_14_0 as an IC SNe; the upper limit assumes that the IC stars are universally old) with the results of Gonzalez et al. 
(2007), who measured the ICL fraction via direct surface brig htness observation s for a sample of 23 low redshift galaxy clusters (points 
with error bars). Also plotted are the intragroup SN results of McGcc & Baloeh (2010), using both their raw numbers (dashed box) and 
those made after assuming that the IC stellar population is composed of only old stars (dotte d box). The combined SN results confirm the 
declining IC stellar mass fraction as a function of cluster mass seen by [Gonzalez ct al. ( 2 0071 '). 



20 Sand et al. 




R/R20O 



Fig. 8. — The distribution of cluster SN la in the MENeaCS survey, with the IC SNe marked by the dark histogram. The IC SNe are 
more centrally concentrated than the hosted cluster SN la population. 
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TABLE 1 

MENeaCS intracluster 
supernovae 



MENeaCS ID 


IAUC ID 


UT Date 


Telescope/ 
Instrument 


Type 


ZSNID 


Abelll650.9.13.0 
Abell2495-5-13.0 
Abell85_6_08_0 
Abell399.3.14_0 


SN2009hc 
SN2008ih 


2009-12-21.45 
2009-06-18.58 
2009-07-04.60 
2008-11-28.49 


MMT/Hecto 
Gemini/GMOS 
Gemini/GMOS 
Gemini/GMOS 


la-norm 
la-norm 
Ia-91bg 
la-norm 


0.0836 (0.0048) 
0.0796 (0.0032) 
0.0617 (0.0007) 
0.0613 (0.0025) 
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TABLE 2 

Cluster Deep CFHT Stack Properties 



Cluster z c i us Exposure Time M g limit M r limit /(< L m i„) /(< L m i n ) 
(s) (mag) (mag) a = —1.03 a = —1.5 



Abcll 1650 


0.0838 


2280 


-12.47 












2240 




-13.04 


0.0012 


0.0172 


Abell 2495 


0.0775 


1920 


-11.72 


-12.37 


0.0007 


0.0127 


Abell 85 


0.0578 


1920 


-11.15 


-11.68 


0.0004 


0.0091 


Abell 399 


0.0718 


720 


-12.54 












2400 




-12.56 


0.0009 


0.0138 



Intracluster Supernovae MENeaCS 



23 



TABLE 3 

Intracluster SN bias correction factors from § 14.11 



Scenario 


Correction Factor 


Hosted to hostless SN detection efficiency 


0.91 


Spectroscopic availability bias 


0.91 


Hostless follow-up bias 


strike Abell2495_5_13_0 from IC calculations 



24 Sand et al. 



TABLE 4 

Intracluster stellar mass fraction 



Scenario Num. of IC SNe a Num. of hosted SNe a IC stellar mass fraction 



All hosts 

All hosts (with Abell399_3_14_0 as IC) 
Red sequence hosts 

Red sequence hosts (with Abell399.3.14_0) 
Hosts within R < -R200 

Hosts within R < R, 2 oo (with Abell399_3_14_0) 
Red sequence hosts within R < R200 

Red sequence hosts within R < R200 (with Abell399_3_14_0) 



2 
3 
2 


20 
19 

8 


o.oBt° :£ 
o.ia±g : j? 

< 0.29 


3 


7 


< 0.37 


2 
3 
2 


14 
13 
6 


n+ 012 

u ' 11 -0.07 

16+ 013 
u ' iD -0.09 

< 0.36 


3 


5 


< 0.47 



a These are the raw number of SNe, after removing Abcll2495_5_13_0 (sec § |4. 11 . The correction factors described in § |4. 1| arc applied before 
the IC stellar mass fraction is calculated. Sec Equation[2] 



